Astronomy <fc Astrophysics manuscript no. 3339MaH 


February 2, 2008 


(DOI: will be inserted by hand later) 





Searching for galaxy clusters using the aperture mass statistics 

in 50 VLT fields* 

M. Hetterscheidt 1 , T. Erben 1 , P. Schneider 1 , R. Maoli 2 ' 4 , L. Van Waerbeke 3 , and Y. Mellier 4 ' 5 

1 Institut fur Astrophysik und Extraterrestrische Forschung (IAEF), Universitat Bonn, Auf dem Hiigel 71, D-53121 
Bonn, Germany 

2 Department of Physics, University "La Sapienza", P.le A. Moro 2, 00185, Rome, Italy 

3 Department of Physics and Astronomy, University of British Columbia, Agricultural Road 6224, Vancouver, V6T 
1Z1, B.C., CANADA 

4 Institut d'Astrophysique de Paris, UMR7095 CNRS, Universite Pierre & Marie Curie, 98 bis boulevard Arago, 
75014 Paris, France 

5 Observatoire de Paris. LERMA. 61, avenue de l'Observatoire, 75014 Paris, France. 
Received 29 April 2005/ Accepted 30 June 2005 

Abstract. Application of the aperture mass (M ap -) statistics provides a weak lensing method for the detection of 
cluster-sized dark matter halos. We present a new aperture filter function and maximise the effectiveness of the 
M a p-statistics to detect cluster-sized halos using analytical models. We then use weak lensing mock catalogues 
generated from ray-tracing through JV-body simulations, to analyse the effect of image treatment on the expected 
number density of halos. Using the M ap -statistics, the aperture radius is typically several arcminutes, hence the 
aperture often lies partly outside a data field, consequently the signal-to-noise ratio of a halo detection decreases. 
We study these border effects analytically and by using mock catalogues. We find that the expected number 
density of halos decreases by a factor of two if the size of a field is comparable to the diameter of the aperture 
used. We finally report on the results of a weak lensing cluster search applying the M ap -statistics to 50 randomly 
selected fields which were observed with FORS1 at the VLT. Altogether the 50 VLT fields cover an area of 0.64 
square degrees. The /-band images were taken under excellent seeing conditions (average seeing w 0'.'6) which 
results in a high number density of galaxies used for the weak lensing analysis (n ps 26arcmin -2 ). In five of the 
VLT fields, we detect a significant M ap -signal which coincides with an overdensity of the light distribution. These 
detections are thus excellent candidates for shear-selected clusters. 

Key words, gravitational lensing - galaxy clusters 
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. F~ i 1. Introduction dictions. The traditional way of obtaining the number den- 

. sity and the main properties of galaxy clusters is using 

Galaxy clusters are the largest collapsed structures m , , , , .... ,., , . ., , , 

J . ° ..„ . „ . direct observable quantities, like luminosity, temperature 

the Universe and formed due to the amplification of pri- ,,, . , , , , n . c . , , , ., 

. . . . , 1 of the mtra-cluster gas and the Ime-of-sight velocity dis- 

mordial density mnomogeneities and subsequent merg- , , , mu _r j ± t n 

........ . persion of the cluster members. 1 he disadvantage of these 

mg processes. Measuring their distribution and structures , . , , , vc , 

. . .. . methods is that they are based on simplified assumptions, 

over a large redshift interval provides crucial mforma- , , , ... , c . , . , , 

° such as hydrostatic equilibrium of the mtra-cluster gas, 

tion about the history and the large-scale structure of . . , . , / ... 

, TT , _ , . . . . , . vinal equilibrium and/or spherical symmetry. 

the Universe. Because their formation and evolution are 

almost exclusively driven by gravity, the number den- Weak gravi tational lensing provides an opportunity to 

sity, main cluster properties and their dependence on measure the ( pro jected) mass distribution without mak- 

redshift can be predicted analytically and by numerical ing a n y of the assumpt ions mentioned above ^Kaiser et alJ 

simulations (|Press fc Schechter||l974HLacey fc Colej|l993t Ischneiderl ITflfli . Furthermore, it is totally inde- 

| Navarro et alj | l996t 11223 1 Jenkins et alj | 200l|) . There are pendent of the baryonic content. Aside from the anal- 

several observational m ethods to test the theoretical pre- ysis of already known mass concentrations like galaxy 

Send offprint requests to: Marco Hetterscheidt, e-mail: clusters, weak lensing techniques can be used to per- 

mhetter@astro.uni-bonn.de form a blind search for hitherto unknown mass concen- 

* based on observations with FORSIOVLT operated by ESO trations, with which it may then be possible to com- 

(programme 63.O-0039A) pile a purely shear-selected cluster sample. Due to high 
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demands on data quality, only in recent years sev- 
eral groups have started to use this method. So far, 
only a few galaxy cluster candidates have been reported 
in the literature llMivaz aki et a l.l 120021: IWittman et alJ 
200 lll200aiSchirmer et al.ll2003l l2004HDahle et al.ll2003t 



Dietrich et all l2004 - These candidates could be iden- 
tified with overdensities of bright galaxies showing the 
presence of 'regular' clusters. Moreover, four of them 
are spectroscopically confirmed. In addition, three shear- 
selected mass concentrations not asso ciated with an op- 
tical counterpart, have been reported |Erben et alJ feoOO: 
lUmetsu fc FutamaseH200bllDahle et al.ll2003|) . Further in- 
vestigations are necessary to confirm or discard these pos- 
sible cases of mass concentrations with an unusually high 
mass-to-light ratio, given that even one of them would 
have profound impact on our understanding of the evolu- 
tion of dark matter halos an d their baryonic content, see 
IVon der Linden et alJ l|2005(l . 

Assuming a random distribution of galaxy orientations 
in the case of no lensing, a coherent alignment of galaxy 
ellipticities could indicate a mass concentration. A quanti- 
tative way to measure this alignment is the so-called aper- 
ture mass (M ap ) statistics l)S chneiderll 1 9 9 6|) . In this paper 
we analyse the ability of the M ap -statistics to detect mas- 
siv e mass concentratio ns. 

iKruse fc Schneider! l|l999l) performed calculations of 
the expected number density of halos usi ng M an with the 
polyno mial filter function introduced by ISchneider et alJ 
( 1998). In this work we use a more effective filter function, 
which ha s alrea dy b een applied to observa tional data by 
ISchirmerl l|2004h and lSchirmer et all l|2004l) . We calibrate 
the filter function to detect a maximum number density of 
cluster-sized dark matter halos assuming an universal den- 
sity profile. We then use this filter function and apply the 
M a p-statistics to simulations and to a data set obtained 
wi th the VLT. 

IWhite et all l|2002h and lHamana et all l(2004l) used 
numerical simulations to determine the expected num- 
ber density of halos and examined completeness and ef- 
ficiency in a weak lensing survey taking into account 
the noise caused by the ellipticity dispersion of back- 
ground ga laxies and the proje c tion e ffects by large-scale 
structure. Irlennawi k, Spergell l|2005|) introduced the so- 
called 'Tomographic Matched Filtering' scheme which 
combines tomography using redshift information of back- 
ground galaxies and matched filtering. In their work it is 
shown that with photometric redshift information at hand 
it is possible to enhance the number density of clusters 
with high signal - to-noi se ratio significantly. In the paper 
bv lMaturi etaTl l)2004l) a nice derivation of a filter func- 
tion is given with the aim of separating dark matter ha- 
los from spurious peaks in weak lensing maps caused by 
large-scale structure lensing. They use numerical simula- 
tions to show qualitatively the sensitivity and reliability 
of this new filter function and compare the results with 
the conventionally used polynomial filter function. 

In the present paper we describe the creation of syn- 
thetic images from ray-tracing through A-body simula- 



tions. In addition to previous work we use these synthetic 
images to analyse how border effects, image treatment 
(galaxy detection and their ellipticity determination, PSF 
correction and shear estimation) and weighting affect the 
signal-to-noise ratio of peaks and their expected number 
density in the weak lensing maps. We finally report on the 
results of a cluster search applying the M ap -statistics to 50 
VLT fields and compare these to the simulations. Criteria 
are presented with which peaks resulting from real clus- 
ters can possibly be distinguished from noise peaks in the 
weak lensing maps. 

The paper is organised as follows. In Sect. |21 we de- 
scribe the Map-statistics and the new filter function, and 
calculate the detectability and number density of halos. 
Numerical simulations are used in Sect.|2|to create syn- 
thetic images in order to determine the expected number 
density of halos and to study the effect of image treatment 
on the signal-to-noise ratio of peaks in the weak lensing 
maps and the resulting change in the expected number 
density of halos. In Sect. 0] we present the creation of 300 
VLT-sized synthetic images with which the border effects 
are studied. The observed VLT fields are analysed in Sect. 
[5] and the number density obtained is compared with ex- 
pectations determined in Sect.0] In Sect-EJall cluster can- 
didates are presented and studied in detail. A summary 
and conclusions are given in Sect.0 Appendix El provides 
the reader with all the formulas we used to calculate the 
number density of halos. 

2. Using the aperture mass for cluster detection 

In the following we use standard lensing notation. For 
a broader introduction to th e topic, see for example 
iBartelmann fc Schneiderl <l200lf) . 



2.1. Introduction 

The gravitational field of a cluster-sized mass concentra- 
tion causes a distortion of the background galaxy images, 
which is revealed as a tangential alignment with respect 
to the centre of the mass concentrat ion. We us e the a per- 
ture mass statistics, introduced bv ISchneiderl l|l996|) . to 
quantify the detectability of cluster-sized dark matter ha- 
los. The aperture mass, A/ ap , is defined as the spatially 
filtered projected mass distribution, k, inside a circular 
aperture of angular radius 9$ at a position 



M ap (0= / d 2 M0)[/(|0-£|), 



(1) 



where U is a radially symmetric continuous weight func- 
tion. Using a compensated filter function of radius 6q, 



d9 6U(9) = 0, 



(2) 



one can express M ap in terms of the tangential shear "f t 
M ap (0 = J d 2 7t (0;OQ(|0-£|). (3) 
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The observable quantity 7t(0;£) = -Re[7(0)e _2i *] is the 
tangential component of the shear at a position 9 — £ = 
\9 — £|(cos<^>, sin^), where 4> is the polar angle of 9 — 
The filter functions Q and U are related through 



dtf' #'£/(#') - 17(0). 



(4) * 



ISchneiderl l)l996(l showed that the variance a c of M ap is 
computable analytically as 



d<9 6>Q 2 ((9), 



(5) 



where cr e is the ellipticity dispersion of galaxies and n is 
the number density of galaxies in the considered aperture. 
The signal-to-noise ratio (snr) is then 



AL 



snr 



a p 



— / d 2 8 lt (8)Q(9) 
J°° dd9Q 2 (9) 



7T(t; 



(G) 



The main advantages of using the aperture mass statistics 
to search for clusters are that M ap can be derived from the 
shear in a finite region and M ap is not influenced by the 
mass-sheet degeneracy (both points follow from the fact 
that the filter function U is compensated) . In addition, the 
application of M ap to observational data is straightforward 
(one has to change the integral in Eq. Q into a sum over 
galaxy images) and the error analysis is simple, see Eq. 



2.2. An adapted filter 

In ISchneider et alJ ( 1998) a family of polynomial filter 
functions which fulfil the conditions J2Jl and (0J were in- 
troduced for mathematical convenience, 



U(0;0 ) 



q + 2) 2 

7T0 2 



1 + 2 



(0/0v? , (7) 



which corresponds to 



Q{0;6 ) 



(1 + 0(2 + 



7T0 2 , 



(6/6 ) 2 (1 - (0/9 o ) 



(8) 



where 8 is the projected angular distance on the sky from 
the aperture centre, 9q is the filter radius. Throughout the 
paper we will choose 1 = 1. However, in order to find the 
maximum number of dark matter halos one should use 
a filter function Q which maximizes the snr. According 
to the Cauchy-Schwarz inequality, the optimal choice for 
radial ly symmetric halos is Q(9) cx 7t(#), see ISchneiderl 
lll99fih . 

Assuming the univers al density pro fi le found b y 
Navarro, Frenk & White ^Navarro et all ljl99fil Il997h . 
NFW-profile in the following), a reasonable choice for the 
filter functio n Q for observa tional purposes has been in- 
troduced bv lSchirmerl l)2004|) . 



Q(x) = (1 + c a ~ bx + 



-c-\-dx\ 



tanh(cc/a: c ) 
ir9$(x/x c ) 



(9) 




Fig. 1. Left panel: the polynomial filter. Right panel: halo- 
filter for different parameters x c with exponential cut off at 
both ends. The quantity x is the normalized filter radius. 
The maximum of all U- filters is normalized to 1. 



with x := 0/9q. The filter function, which we refer to 
as 'halo-filter' in the following, approximately follows the 
tangential shear profile of an NFW-halo over a large ar- 
range and is mathematically simple. We choose the values 
a = 6 and b = 150 so that Q exponentially drops to zero 
at x = 0. Furthermore, a choice of e.g. c = 47 and d = 50 
results in an exponential cut-off around x = 1. The filter 
function has the nice properties that it downweights the 
inner part of a cluster which is often associated with bright 
galaxies (so no faint background galaxies are visible in the 
centre) and that smooth weak lensing maps are obtained. 

The parameter x c changes the shape of the filter in 
such a way that more weight is placed at smaller radii for 
smaller values of x c . In the right panel of Fig. ^ the filter 
function Q and the corresponding filter U are shown for 
different values of x c . In the left panel o f this figure the 
polyno mial filter functions introduced bv ISchneider et, all 
( 1998J) are displayed for comparison. 

In order to obtain predictions for the expected signal- 
to-noise ratios for halos of mass M at redshift z\, we have 
to specify a mass model fc(i?, Z\, Zs) and a distribution for 
the source redshifts p(z s ). The expected Af ap -signal is then 
given by 



M ap = 



d 2 $ 



dz s K,(-d,z h z s )U(-d)p(z s 



(10) 



The noise is simply given by a c in JSJ. To calcu- 
late the mass profi le we closely follow the work of 
lHamana et alJ ll2004h. We consid er a truncated NFW- 
profile, see Tak ada fc Jainl <|2003), and use the con cen- 
tration parameter introduced by Bullock et al 1 l|200lh . In 
Appendix A all equations to calculate the surface mass 
density n are listed. 

To use the filter function ^ we calculate how the pa- 
rameter x c changes the snr for different halo masses, red- 
shifts and filter radii; see Fig.|3for an example. It turns out 
that a good choice for all reasonable combinations of halo 
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Fig. 2. The dependence of the snr on the filter parame- 
ter x c for different filter radii 9q for two different cluster 
masses and redshifts. Left panel: high cluster mass, low 
redshift. Right panel: low cluster mass, high redshift. 



masses, redshifts and filter radii is x c — 0.15. Fixing x c , 
and using the filter radius 9q = 6' (this filter radius is the 
optimal choice to find the maximum number of halos in 
the redshift range z S [0; 0.95], see Fig.HjJ we calculate the 
snr for different masses and redshifts of halos, given the 
redshift distribution o f background galaxies introduced by 
iBrainerd et"all ljl99fih . 



p(z) 



T[{l + a)/[3]z \z 



— exp 



(11) 



where T 
and zn 



is the Gamma-function, and a 



2, (3 



0.8, which results in a mean redshift (z) 



1.5 

1-2 

/llson c t a"fll200lh . Furthermore, we assume in the fol- 
lowing a typical weak lensing survey with a number den- 
sity of galaxies of n = 30 arcmin -2 and an ellipticity dis- 
persion of <7 e = 0.4. The results are shown in Fig. [31 With 
the halo-filter function we should be able to detect halos 
with a signal-to-noise ratio larger than 4 down to masses 
of 1O 14 M for redshifts lower than 0.3. 



2.3. Number density of halos 

To calculate the number density of significant peaks in 
the aperture mass map above a given snr-threshold, N(> 
s nrt), resulting from real ha los, we closely follow the work 
of lKruse k, Schneider! l)l999(l . It is assumed that dark mat- 
ter halos are distributed according to the Press-Schechter 
theory . We use the fitting formulae given in lNavarro et alJ 
( 1997) to compute the number density N of halos. The 
aperture mass is a monotonically increasing function of 
halo mass M for constant values of the lens redshift z\ 
and filter radius 9q, therefore it can be inverted for a given 
threshold value M* p . As the noise in Eq. (JBJ only depends 
on the filter function in use, the snr-threshold value is 
snrt — Mlp/cc- The number of halos in a given proper 




Fig. 3. Detectability of mass concentrations with an 
NFW-profile for different masses M and redshifts z\. The 
contour lines indicate the snr. We assume J7a = 0.7, 
f2 m = 0.3, number density of galaxies n — 30/arcmin 2 
and an ellipticity dispersion of <t c = 0.4. The filter radius 
is chosen to be do = 6' and x c = 0.15. 

volume with a mass greater than Mt(snr t , z\, 9q), and thus 
a snr greater than snr t , is given by, 



N(> snr t ) = / dVp(l + zi) 



dMjV halo (M, *l) (12) 

M t 

l|Kruse fc SchnfiidfiT<ll99d> . The quantity iV halo dKdiW is 
the number of halos in the comoving volume dV c with 
mass in the interval dM. We now calculate the number 
density 7V(> snrt) of halos above a given snr-threshold, 
assuming the mass profile and observational parameters 
(n = 30/arcmin 2 , <j e = 0.4, (z) = 1.2) of Sect. l2~2l 
As can be seen in Fig. the theoretical number den- 
sity A^(> snr t = 4) calculated for the halo-filter is w 3 
ti mes larger than for t h e po lynomial filter. In comparison 
to iKruse fc Schneider! l|l999» our results for N{> snrt) 
using the polynomial filter are much lower, because our 
assumed ellipticity dispersion is twice as large, our mean 
redshift is lower and we use a slightly different mass pro- 
file. 

Taking into account only the noise induced by the in- 
trinsic ellipticity distribution of the background galaxies, 
we assume that the difference between the real value of 
M ap and the measured M ap follows a Gaussian distribu- 
tion, 



p(AM ap ;0 o) 



exp 



AM 2 p 
2 a 2 



(13) 



The observable number density N of halos above a given 
snr-threshold is then obtained by convolving the theoreti- 
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Table 1. The number of halos per 10 square degrees with 
a signal-to-noise ratio larger than 4 and 5 for the theoret- 
ical number density, N, the observable number density, N 
(both calculated with a filter radius of 6 = 6') and the 
observable number density, for which it is assumed that 
the area is analysed with different filter radii, iV var i. 



2i G [0,0.95] 


iV/10deg 2 


TV/10 deg 2 


7V vari /10deg 2 


snr > 4 


37 


55 


61 


snr > 5 


16 


23 


25 



Table 2. The maximum number of halos with a signal-to- 
noise ratio larger than 4 per 10 square degrees per redshift 
interval. The value 8q is the filter radius which maximises 
the number density for the given redshift interval. 





A> max [lOdeg" 2 ] 


8o [arcmin] 


[0.0,0.15] 


4.6 


13 


[0.15,0.25] 


14.7 


8 


[0.25, 0.35] 


14.7 


6 


[0.35, 0.45] 


12.3 


6 


[0.45, 0.95] 


14.4 


5 



cal nu mber density N by p(AM ap ), see lKruse fc Schneider! 
()l999(l . This results in an increase in the number den- 
sity of peaks of high significance, see Fig. 21 and Table 
n This, however, is only an approximation. The noise 
not only changes the peak height but also the peak po- 
sition. A neighbouring pixel of the original peak maxi- 
mum (without noise) can be higher after adding noise, so 
the given calculations are just an upper limit. A profound 
analysis of the influence of Ga ussian random fields o n the 
peak statistics can be found in IVan Waerbekd l)200Cf) and 
Jain fc van Waerbekd l|2000|) (however, for high snr our 
approximation is sufficient). 

The values of projected angular radii on the sky of 
galaxy clusters having various virial radii and redshifts 
are very different. Considering this fact one applies, in 
practice, the aperture mass statistics with varying filter 
radii. To estimate the increase of the number density by 
applying different filter radii, we determine, for different 
redshift bins, the aperture radius for which we obtain the 
maximum number density of halos (Fig. [BJ) . The maximum 
number is then added up (Table ^ 13 • We find that the 
expected number density of halos with a signal-to-noise 
ratio larger than four exceeds 61 per 10 square degree, 
which is only slightly larger than using a fixed filter ra- 
dius of 6*o = 6' (55 per 10 square degree). This is due to 
the fact that the radius for which one expects the maxi- 
mal number of halos is approximately the same for every 
redshift interval (Fig. • 

In the real world, not every galaxy cluster is relaxed 
and has an NFW-profile. It is therefore difficult to esti- 
mate the difference in the number density between a fixed 




3 4 5 6 



snr 

Fig. 4. Theoretical number, N, and the observable num- 
ber, N, of halos per square degree with a signal-to- 
noise ratio larger than snr, calculated for the halo-filter 
(#o = 6') and, as a comparison, for the polynomial fil- 
ter (9o = 4'). The observable number density is obtained 
by convolving the theoretical number density with the 
Gaussian distribution p(AM ap ; 9 ) (|13|) . 

and variable filter radius because our filter function is opti- 
mised for NFW-profiles. This is worth examining in detail 
with numerical simulations in a future paper. 

In Table [21 it can also be seen that even for redshifts 
larger than 0.45 we expect to find about 15 NFW-halos 
with a signal-to-noise ratio larger than 4 in a 10 deg 2 sur- 
vey. 

2.4. M ap applied to images 

We now describe the application of the M ap -statistics to 
images. It is straightforward to construct an unbiased es- 
timator M ap for the integral by a discrete sum over 
observed galaxy ellipticities e t ■ Considering the coordinate 
origin to be at the centre of the aperture this then reads 
as 

H n A — ' 

% 

where n is the number density of galaxies in the consid- 
ered aperture and et = — Re[e(i9)e~ 2lc ^] is the tangential 
ellipticity. The discrete dispersion Cd of M ap in the case 
of no lensing can be calculated by squaring Eq. I|14f> and 
taking the expectation value, which leads to 

^ = |E^). ( 15 ) 
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Fig. 5. Dependence of the number of halos per square de- 
gree with a signal-to-noise ratio larger than 4 and larger 
than 5, respectively, on the aperture radius of the halo- 
filter. The optimal choice for the halo filter radius to find 
the maximal number of halos with a snr larger than 4 is 
Oo = 6'. 




6 8 10 
Q Q [arcmin] 



6 8 10 
6 Q [arcmin] 



Fig. 6. Dependence of the number of halos per square de- 
gree with a signal-to-noise ratio larger than 4 on the aper- 
ture radius for different redshift bins. Left panel: theoret- 
ical number density (N). Right panel: observable number 
density (N). 

where cr 2 — (|e| 2 ) is the ellipticity dispersion and n the 
number density of galaxies. We used the fact that the el- 
lipticities of different images are not correlated ((euetj) = 
8ij(r1/2). In the case of weighting, the discrete aperture 
mass M' changes to 



M. 



and the discrete dispersion then reads 

2 .^Eil^MWO 



(16) 



(17) 



2(E^i) 

see ISchirmerl l|2004h . The weighting factors iw, are calcu- 
lated from the uncertainty of e t [see Eq. (154*1) ]. A regular 



grid of aperture centres is now placed over a data field 
and the snr is calculated for every grid point. In this way 
weak lensing maps (M ap -maps) are obtained for every im- 
age in which mass concentrations are revealed as peaks in 
the map. 



3. Applying M ap to numerical simulations 

In this section we investigate the ability of the aperture 
mass statistics to detect mass concentrations by applying 
M ap to images obtained from ray-tracing through ACDM 
N-body simulations (£!a = 0.7, f^o = 0.3, a s = 0.9, h = 
0.7). These simulations were kindly made availab le by 
Takashi Hamana [details see lHamana et aP l|2004l) ]. We 
create twelve initial catalogues of randomly distributed 
galaxies using the programme stuff 1 (E. Bertin). For a 
detailed description of the galaxy morph ology and magni- 
tude distribution see lErben et al.l (2001). The galaxies are 
assumed to be at a fixed redshift z = 1 and are sheared 
according to the shear map of the ray-tracing simulations, 
meaning that we modify the intrinsic galaxy ellipticity e 
by the shear 7 present at that position. In the following 
these galaxy catalogues are called input catalogues. 

The input catalogues are used to create synthetic im- 
ages using the programme SkyMaker 2 by E. B ertin. A 
short description is given in lErben et al ] ll200lh . Twelve 
30' x 30' images resulting in a 3 deg 2 survey are obtained 
for the twelve catalogues. These images are treated in ex- 
actly the same way as real data (like object detection, 
PSF correction, same cuts, weighting), see Sect. 15.21 For 
the object detection we u tilise the programme SExtractor 
l|Bertin fc Arnoutsl Il996[) . An object is detected if 3 
contiguous pixels are la above the sky background 
(SExtractor parameter settings DETECT_MINAREA=3, 
DETECT_THRESH=1). The obtained galaxy catalogues 
are called output catalogues. 

To exclude the effect of false detections we only take 
into account those objects which are present in both cat- 
alogues. The mean ellipticity dispersion of the galaxies of 
the input and output catalogues is tr £ = 0.32 and the mean 
galaxy number density is n — 19arcmin~ 2 . Note that the 
galaxy number density quoted is the number density of 
background galaxies, as all galaxies are placed at redshift 
z = l. 

In the following we apply the M ap -statistics using the 
halo-filter for two different filter radii (8q = 1000 pixels = 
3'.8 and #0 = 1500 pixels — 5'.7) on the input and output 
catalogues as described in Sect. 12.41 where the mesh size 
of the grid is 10" x 10". We examine various effects on 
the snr of peaks and on the number density of signifi- 
cant peaks (snr > 3) in the M ap -maps. At this stage it is 
worth mentioning that significant peaks in M ap -maps can 
originate from 



Available at: 
http://terapix.iap.fr/cplt/oldSite/soft/stuff/ 

2 Available at: 
http://terapix.iap.fr/cplt/oldSite/soft/skymaker/ 
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1. massive halos, 

2. low mass halos which are boosted to high snr by noise, 

3. pure noise peaks, 

4. projection effects (two or more halos at different red- 
shift along almost the same line of sight are merged in 
the M a p-map due to the smoothing effect of the filter 
function) , 

5. real substructure in massive halos, 

6. false high peaks in the vicinity of a massive cluster due 
to noise. 



In this work we cannot distinguish between these cases. 
Nevertheless, we can estimate the number of pure noise 
peaks by randomising the background galaxies in each 
field, see Sect. 13.41 The influence of points 4-6 can be re- 
duced if we consider only peaks with high signal-to-noise 
ratios (snr > 3). 

For the analysis of various effects on the snr of peaks 
in the M ap -maps we investigate the maps by eye. For the 
statistics of the number density of halos we define a peak 
in the M ap -map as a pixel with maximum value compared 
to the surrounding 24 neighbour pixels. By applying this 
definition and a mesh size of 10" x 10" we reduce the 
number of false peaks in the vicinity of a massive cluster 
with snr > 3 (point 6) . The difference between the number 
density of peaks found according to this definition and the 
investigation by eye is insignificant. 



3.1. Example of the simulation 

In Fig. we present one of the re-maps from the N-body 
simulations and compare it to the snr-maps for different 
filter functions and filter radii. The most prominent fea- 
tures in the re-map (upper left) are detected with a high 
significance in the sror-maps independent of the filter in 
use. The halo-filter not only detects mass concentrations 
with higher significance compared to the polynomial filter, 
but also resolves the maxima in the re-map, indicated by 
the white circles in Fig. 7. 



3.2. Two different halo-filter radii 

We apply the M ap -statistics for two different halo-filter 
radii (f?o = 3'.8 and 6*o = 5'.7) to the input and output cat- 
alogues of the numerical simulations. As expected from the 
calculations in the previous sections, the snr on average 
increases for larger radius, see lower left diagram in Fig. 
HT)I as example for the output catalogue. With the larger 
radius more significant peaks are detected for the input 
catalogue. This is illustrated in the upper right diagram 
of Fig. ^3 11 peaks per square degree are detected with a 
snr larger than four using the halo-filter with a radius of 
9q = 5'.7, compared to only 6.5 using the halo-filter with 
a radius of 9q = 3'.8 (see also Tab.EJ- 






Fig. 7. An example of the N-body simulations. Upper 
left: re-map, upper right: snr-map, polynomial filter 
with 9q = 3'.8, lower left: snr-map, halo-filter with 
9q = 5'.7, lower right: snr-m&p, halo-filter with 9o = 5'.7 
after source extraction and PSF correction. Independent 
of the filter in use the most prominent features in the Ac- 
map are detected with high significance. White circle: in 
contrast to the polynomial filter the halo-filter resolves the 
three peaks in the re-map. The size of each field is 30' x 30'. 
Contour lines in kappa-map: 0.05, 0.1, 0.15 and 0.2; in all 
snr-maps: 2, 3, 4, 5, 6. 

Table 3. Number of peaks per one square degree with a 
signal-to-noise ratio larger than four. 





N(snr > 4) 


halo-filter, input cat., 9q = 5'.7 


11 


halo-filter, output cat., 9o = 5.7 


6 


halo-filter, output cat., weighting, 8o = 5'. 7 


7.5 


halo-filter, input cat., border effects, 9o = 5'. 7 


6 


halo-filter, output cat., border effects, 9o = 5'. 7 


3 


halo-filter, input cat. 9q = 3'.8 


6.5 


halo-filter, input cat., border effects, 9q — 3'.8 


4 


halo-filter, output cat., border effects, 9o = 3'.8 


3 


poly-filter, input cat., 8o = 3'. 8 


2.5 


VLT images, 9 = 5'.7 


2 


VLT images, 9o = 3'.8 


3 



3.3. Comparison between halo- & polynomial filter and 
the analytical model 

The snr of peaks in the M ap -maps obtained from the in- 
put catalogues is significantly lower if the polynomial fil- 
ter is used, see lower right diagram in Fig. As a con- 
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sequence, the number density of peaks in the M ap -maps 
using the polynomial filter for a given snr is significantly 
lower compared to the halo-filter. We obtain 11 peaks per 
square degree with a snr larger than four for the halo- 
filter (6> = 5'.7), but only 2.5 for the polynomial filter 
(8q = 3'.8), see lower right diagram in Fig. 1111 

In Fig. 1111 the number densities of dark matter ha- 
los for the halo-filter and the polynomial filter obtained 
from the analytical model are also plotted. In this sec- 
tion the same parameters are used for the calculations as 
for the simulations (fixed redshift z = 1, 6q = 5'7 for 
halo-filter, 6q = 3'.8 for polynomial filter, cr e = 0.32 and 
n = 19/arcmin 2 ). Notable is the difference between the 
analytical model and the simulations. For snr less than 
ps 4 the reasons are that peaks in the synthetic data can 
also originate from the effects listed in the introduction 
of Sect. |3 The main reason for the difference for high 
snr is that more massive clusters are present in the syn- 
thetic data (note that the area of the data only covers 
3 deg 2 ) as is indicated by the jump in the number density 
at snr = 5.6 (halo-filter curve) which is then carried over 
to lower snr. Another reason could be that we have used 
the Press-Schechter model to calculate the number den- 
sity oflialos J _2diichj ; inderpredicts the more massive ones, 
see l.Tenkins et alJ i|200lh . 

From the calculations we expect to detect 3 times more 
halos with the halo-filter than with the polynomial filter 
for a snr > 4. But as the halo-filter resolves substructure 
(Fig. due to the fact that the halo-filter is much nar- 
rower than the polynomial filter) the number density of 
Map-peaks obtained with the halo-filter is even 4.5 times 
larger than the number density obtained with the polyno- 
mial filter. 

3.4. Noise peaks in weak leasing maps 

The number density of halos is contaminated by noise 
peaks caused by a chance alignment of background galax- 
ies. To quantify the number density of these noise peaks, 
we randomise the orientation of background galaxies, ap- 
ply M ap and repeat this procedure 20 times for each of 
the 12 output catalogues. This is done for the two differ- 
ent filter and for two different filter scales, see Fig [S] and 
Tab.0] As expected, the number density of noise peaks of 
a given snr-threshold is lower for a larger filter radius and 
is larger for the halo-filter compared to the polynomial 
filter, because it is narrower than the polynomial filter 
and therefore puts a high weight to a smaller number of 
galaxies. 

The number density of pure noise peaks is then sub- 
tracted from the total number density of peaks to estimate 
the number density of 'real' peaks that are due to a real 
overdensity in the K-map of the simulated field (Fig. [HJ) ■ 

We now answer the question of how many real peaks 
in the weak lensing maps we obtain if we only allow a 
contamination ratio (number of noise peaks to total num- 
ber of peaks in the weak lensing maps) of 20%. We plot 



Table 4. Number of noise peaks per one square degree. 





N(> 3) 


iV(>4) 


sim., halo-filter, 8q = 5'.7 




17 


0.75 


sim., halo-filter, 8o = 3'.8 




30 


1.28 


sim., halo-filter, border eff., i 


o = 8.7 


21 


0.81 


sim., poly, filter, 6q — 5'.7 




3.7 


0.10 


sim., poly, filter, 9o = 3'.8 




7.7 


0.22 


sim., poly, filter, border eff., 


% = 8.7 


3.4 


0.10 


VLT images, halo-filter, 8o = 


8.7 


20 


0.68 


VLT images, halo-filter, 8o = 


8.8 


30 


0.85 



this ratio in Fig. If we assume a contamination of 20% 
we obtain 11.6 real peaks for the halo-filter and 8.9 real 
peaks for the polynomial filter. If we take into account the 
contamination ratio, the difference of the efficiency of the 
two filter types is much less pronounced compared to the 
efficiency obtained if we only take into account a given 
snr-threshold, see last section. In other words, to com- 
pare the efficiency of different filters it is not sufficient 
to compare the number density of peaks for a given snr- 
threshold, but to compare the number density for a given 
contamination of noise peaks. 

3.5. The effect of image treatment 

Image treatment includes galaxy detection and the de- 
termination of galaxy quadrupole moments, PSF correc- 
tion (anisotropy- and P g -correction) and catalogue filter- 
ing. This will be described in detail in Sect. 15.21 In this 
work we discuss the effects of image treatment on the snr 
and the change in the number density of M ap -peaks as a 
whole. Detailed insights into the impact of different steps 
of the image treatment or KSB algorithm and its imple- 
mentation on shear estimates will be presented for two 
different KSB-pipclines (Bonn and Edinburgh pipeline) in 
Hetterscheidt et al. (in prep.). 

We apply the M ap -statistics to the input and output 
catalogues and compare the change of the snr of peaks in 
the M a p-maps, see upper left panel in Fig.EH It is clearly 
visible that the image treatment lowers the snr of peaks 
in the M ap -maps of the input catalogue significantly. 

We quantify this effect by plotting the number den- 
sity of snr-peaks in the M ap -maps, see Fig. El The lower 
left panel of this figure compares the number density of 
peaks in the M ap -maps of the input catalogue and the 
output catalogue. The image treatment definitely lowers 
the number density of peaks in the M ap -map, the number 
density of the output catalogue is lower than that of the 
input catalogue by a factor of two. We obtain 11 peaks 
per square degree with a snr larger than four for the in- 
put catalogue, but only 6 for the output catalogue (see also 
Tab. 0). The ratio between the number density of peaks 
obtained by analysing the input catalogue and output cat- 
alogue can therefore be reduced using more conservative 
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Fig. 8. Comparison of number density of peaks in the 
Af a p-maps per square degree with a signal-to-noise ra- 
tio larger than snr of the output catalogue for halo- and 
polynomial filter and different filter radii. Solid lines with 
error bars: total number of peaks per square degree in 
the Map-maps. Dotted lines: number of noise peaks per 
square degree, resulting from randomisation of the orien- 
tation of background galaxies. Bold solid lines: number of 
'real' peaks (difference between the total number of peaks 
in the M ap -maps and number due to randomisation) per 
square degree. The error bars are due to Poisson statistics 
obtained from 3 square degrees. For clarity they are only 
plotted for the total number density. 



SExtractor parameter settings (like a larger number of 
contiguous pixels). 

Faint galaxies have intrinsically the same size as the 
PSF or smaller. Hence, the observed galaxy images are 
composed of only a few pixels and the cllipticity determi- 
nation via quadrupole moments of the surface brightness 
is extremely noisy as is the PSF correction which is also 
calculated by means of the quadrupole moments. The re- 
duction of the snr of peaks and the number density of 
peaks therefore follows from the fact that shear informa- 
tion is destroyed due to the noisiness of galaxy images and 
the subsequent extremely noisy correction process. 

Very long exposure times with ground-based telescopes 
would not increase the number density of background 
galaxies appreciably because additional fainter galaxies 
are much smaller then the PSF so that one cannot mea- 
sure the pre-seeing surface brightness properly. For future 
large weak lensing surveys it is therefore not efficient to 
propose for a few very deep images to perform a weak 
lensing analysis compared to many shallower images. 



Fig. 9. Contamination of noise peaks. Ratio between num- 
ber of noise peaks and total number of peaks for the halo- 
and polynomial filter and for two different filter radii. The 
number of real peaks is the difference of the total number 
of peaks and the noise peaks at the snr of 20% noise peak 
contamination. 

3.6. The effect of weighting 

To study the effect of weighting (the weighting scheme is 
explained in Sect. I5.2")l we apply the aperture mass statis- 
tics to the output catalogues without weighting, Eq. l(T3jl. 
and with weighting, Eq. I|16|) . It turns out that weighting 
in the form described in Sect. l5.2l on average raises the snr 
of peaks in the M ap -map only slightly, see upper right dia- 
gram in Fig. llUI The small influence of weighting is mainly 
due to the fact that we exclude PSF corrected ellipticities 
of more than 0.8. We plot the number density of peaks 
in the M ap -maps with and without weighting in the up- 
per left diagram of Fig. 1111 The number density increases 
slightly using weighting. We obtain 6 peaks per square de- 
gree with a snr larger than four for the output catalogue 
without weighting, but 7.5 for the same catalogue with 
weighting. 

4. Border effects 

A maximum snr for a given filter function Q is attained 
when the centre of a radially symmetric mass concentra- 
tion and the circular aperture centre coincide. However, if 
the aperture lies partly outside a field, the snr is lowered 
due to the fact that the number of galaxies from which 
the shear is estimated decreases. To avoid these border ef- 
fects, the distance x of the centre of the circular aperture 
from the edge of the field must be larger than the aper- 
ture radius. The resulting subfield would be significantly 
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snr of output (0 O =3.8') snr of input (halo; 8 =5.7') 

Fig. 10. Comparison of the signal-to-noise ratio of peaks 
in the Af ap -maps for different analyses, upper left: 
Comparison between input catalogues and output cata- 
logues (after image treatment), upper right: Comparison 
between weighting and no weighting (for output cata- 
logues), lower left: Comparison between two different 
halo-filter radii, lower right: Comparison between halo- 
filter and polynomial filter. 



smaller, even for wide-field images (the radius of the halo- 
filter, for which one obtains the maximum number density 
of NFW-halos is 6 « 6', the subfield of a 30' x 30' field 
would then be 18' x 18', consequently 64% smaller). 

At this point it is worth mentioning that the condition 
of a compensated filter function U [see Eq. (|2J|] is no longer 
valid if the circular aperture lies partly outside the field, 
so that the M ap -value measured is not fully related to the 
projected mass distribution, k. Note, however, that the 
goal of this method is to find mass concentrations and not 
to get detailed information about the mass distribution. 




snr snr 



Fig. 11. Number density of peaks in the M ap -maps per 
square degree with a signal-to-noise ratio larger than 
snr obtained from ray-tracing through N-body simula- 
tions. Upper left: dotted histogram: number density 
of peaks in the M ap -maps of output catalogue (9 = 5'.7; 
no weighting). Solid histogram: same but with weighting. 
Upper right: solid histogram: number density of peaks 
in the M ap -maps of input catalogue (9$ = 5'7; no weight- 
ing). Dotted histogram: same but with a filter radius of 
0o = 3'.8. Lower left: solid histogram: number density 
of peaks in the M ap -maps of input catalogue (9a = 5'.7; 
no weighting). Dotted histogram: same but for the output 
catalogue. Lower right: solid histogram: number density 
of peaks in the Map-maps of input catalogue for the halo- 
filter (9o = 5'.7; no weighting). Dotted histogram: same 
but for the polynomial filter with 9o = 3'.8. Long dashed 
line: analytical model with same parameters as simula- 
tions (halo-filter) . Long-short dashed line: same but for 
the polynomial filter. The error bars are due to Poisson 
statistics obtained from 3 square degrees. Because of clar- 
ity they are not plotted for every line. 



4.1. Border effects for an SIS profile 

As an example we calculate how the snr of a cluster with 
a profile of a singular isothermal sphere (SIS) centred in 
the circular aperture changes as a function of the distance 
x of the aperture centre to the edge of the field (the cal- 
culations below are done with the polynomial filter func- 
tion, because in this case the solution is analytic). We 
introduce the parameterised radius #(</?) = xj cos tp for 
<p € [—tpi,tpi] with ipi = arccos(x/#), see Fig. El With 
this we then obtain for the signal 



where 



I Ml 



and 



IM2 — 



dip 



pi 



dip / dl?'T?'7t(0')W)- 



(19) 



(20) 



We solve the integral for an SIS-profile and the filter func- 
tion JSJ|. The tangential shear profile of an SIS reads 



M ap = 2(/a/i + Imi) 



(18) 7t 



0E 

29 ' 



(21) 
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Fig. 12. The figure displays a part of the circular aper- 
ture with radius 0, which lies partly outside a field (shaded 
area). The parameter x denotes the distance between the 
aperture centre and the edge of the field, is the pa- 
rameterised radius and (pi is the integration limit. 



0.6 - 



0.5 




n_=0.3, S)„=0 

halo-filter, M=lxl0"M ai 

halo-filter, M = 2xl0 14 M a , 

poly-filter, M=lxl0 14 M„ 

poly-filter, M 2x10" M„ 

poly-filter, SIS 



where 0e is the Einstein radius. We obtain, 



-0.4 

Ve 



Imi = 



with 



4Ott0 6 



9 2 (U9 2 - 9x 2 



f ± = (±20 9 2 x 2 T 9x 4 ) In <^ cos 



1 - 



2 V1 



X" 



± sin ( -ipi 



and 



Im2 = 



29 E/ 



+ f_ + /+ > , (22) Fig. 13. Displayed is the relative change of the snr if a 
J circular aperture lies partly outside a field and a cluster 

with NFW-profile is positioned at the centre of the aper- 
ture for different filter types, halo masses and halo red- 
shifts. The same is displayed for an SIS-profile, where the 
polynomial filter is used. The quantity x is the distance 
of the aperture centre to the edge and 9 is the radius of 
the aperture. Negative x: aperture centre inside the field, 
positive x: aperture centre outside the field. 



(23) 



(24) 



For the M ap signal of the whole aperture we obtain, M ap = 
40 E /(5 0). For the noise we find, 

2 

<a P = 2^(/ CTl +/ CT2 ), (25) 
where 

I*i = dtp dW'Q 2 (26) 

Jo Jo 

and 

h2 = I dip I dddQ 2 . (27) 

J <pi Jo 

The noise is independent of the signal. We obtain 

^^u^vf^' (28) 

with 

p = 55 9 s + 90 9 6 x 2 + 296 9 4 x 4 - 592 9 2 x 6 + 256 x 8 (29) 
and 

3 



5tt 2 (9 2 



(ic-tpi). 



(30) 



For the snr and its dependence on the distance x from the 
edge of the field, we only have to calculate M ap /<jMa P - The 



result is shown in Fig. 1131 As expected, the snr is reduced 
by a factor of l/y/2 if the centre of the aperture lies at 
the edge. Having only half of the aperture is equivalent to 
having the half of the galaxies in the whole aperture. We 
also calculate, for an SIS-profile and the halo-filter, how 
the snr changes if the circular aperture reaches the corner 
of a field. The result is shown in Fig. 

4.2. Border effects for an NFW profile 

We calculate numerically how these border effects affect 
the snr of a halo with NFW-profile positioned in the cen- 
tre of the aperture. The results are also shown in Fig. 
IT3l We conclude that if one uses the halo-filter the snr 
drops significantly (more than 5%) only if the centre of 
the aperture is closer to the border of a field than 0.4 x 9q. 

4.3. Influence of border effects on the number density 
of peaks in the M &p -map 

Using the filter function J!2J), the optimal aperture radius 
to detect a maximum number of NFW-halos is theoret- 
ically O = 6'. As each VLT/FORS1 field only covers 
6'.8 x 6'.8, the optimal radius has the same size as the en- 
tire VLT field. We have seen that the snr decreases if 
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Fig. 14. Relative change of the snr if the circular aperture 
reaches the corner of a field. Assumed is an SIS-profile 
and the halo-filter. The contour levels from lower left to 
upper right are 95% to 5% of the full aperture snr. The 
centre of the cross indicates the corner of the field, where 
Xi/Oq — X2/O0 = 0; $o is the aperture radius. 



using the aperture mass statistics in 50 VLT fields 




3 4 5 3 4 5 



snr snr 

Fig. 15. Influence of the border effects on the number den- 
sity of halos. Displayed is the comparison of the number 
density for twelve 30' x 30' images (solid lines) and for 
300 6' x 6' images (dashed lines) of the input catalogue. 
Left: halo-filter radius 9q = 5'.7. Right: halo-filter radius 
(?o = 3'.8. For a better comparison the error bars are only 
plotted once. 



The number of real peaks is therefore significantly smaller. 



the aperture is partly outside a field. We now determine 
how this affects the expected number density of halos. 
Therefore, we split each of the 12 fields of the numerical 
simulations into 25 subfields all having a size of 6' x 6', 
obtaining 300 VLT-field-sized images. We apply the Map- 
statistics to these subfields for two filter radii (60 = 3'.8 
and 9q = 5'.7) and find that the number density of peaks 
drops by a factor of ~ two compared to the large fields, 
see Fig. ED and Table El 

4.4. Influence of border effects on the number density 
of noise peaks 

In this section we study the influence of the border effects 
on the number density of noise peaks and the resulting 
number density of real peaks in weak lensing maps of small 
size. Noise peaks and real peaks are defined as in section 
13.41 In Fig.^Jwe compare the number density of all peaks, 
noise peaks and the resulting real peaks, for twelve 30' x 30' 
images (total fields) and for 300 6' x 6' images (subfields) 
of the output catalogue for two different filter (filter radius 
in all cases #0 = 5'.7). In Tab. 21 two examples are shown. 

The tremendous influence of the border effects on the 
number density of real peaks (difference between total 
number and number of noise peaks) can be seen in Fig. 1161 
The total number of peaks in the subfields drops and the 
number of noise peaks rises compared to the total fields. 



5. Weak lensing analysis of the VLT data 

5.1. The data 

For the current work we observed 50 uncorrelated fields 
with the FORS1 camera at the VLT (UT1, ANTU) at 
Paranal. These fiel ds have been used before for a cos- 
mic shear analysis l(Maoli et al.ll200lh . The observations 
were taken under optimal conditions (average seeing ~ 
Of! 63) in /-band during the Period 63 (March 1999 to 
September 1999). The FORS1 camera is equipped with 
a 2048 x 2048 pixel CCD chip with a pixel size of 0.24 /urn 
which corresponds to ~ 0'.'2 in the standard mode of 
the instrument. Consequently, each field covers 6'.8 x 6'.8. 
Altogether the 50 fields cover ~ 0.64 deg 2 . The total ob- 
serving time for each field was 36 minutes, which resulted 
in a limiting magnitude of Jab ~ 24.5, corresponding to 
a mean galaxy number density of 26 arcmin -2 in the fi- 
nal weak lensing catalogue. The expected mean redshift 
of the lensed sources is (z) w 1 . The fields represent 
50 independent lines-of-sight chosen in such a way that 
they are neither biased towards overdense or underdense 
regions, nor bright stars or other very luminous objects 
are within or close to the field, in order to have enough 
background galaxies for the weak lensing analysis. For the 
following analysis the VLT data was already reduced. For 
further information about the observations, se lection crite- 
ria of the fields and data reduction we refer to lMaoli et alJ 
fcOOll) . 
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Fig. 16. Influence of the border effects on the number den- 
sity of all peaks and noise peaks. Displayed is the compari- 
son of the number density for twelve 30' x 30' images (total 
fields) and for 300 6' x 6' images (subfields) of the output 
catalogue for two different filters (filter radius in all cases 
9q = 5'.7). Solid lines with error bars: total number density 
of peaks in the Map-maps. Dotted lines: number density 
of noise peaks, resulting from randomisation of the orien- 
tation of background galaxies. Bold solid lines: difference 
between total number density of peaks and noise peaks. 

5.2. Catalogue creation and ellipticity correction 

In this section we briefly elucidate our catalogue creation, 
catalogue filtering, PSF correction and weighting scheme. 
Raw catalogue. SExtractor is used to create two pri- 
mary catalogues of all objects in the /-band image which 
consist of at least N = 3 and N — 5 contiguous pix- 
els (SExtractor parameter 'DETECTJV1INAREA') with 
a flux greater than the k — la and k — 2a sky level noise 
(SExtractor parameter 'DETECT.THRESH') for all 50 
VLT fields. 

First catalogue filtering. All objects for which prob- 
lems concerning the determination of shape or position 
occur are rejected (e.g., objects near the border, with neg- 
ative total flux, with negative Qn + Q22 (see below), or 
with negative semi major and/or semi major axis). Two 
catalogues of 107100 objects for N = 3, k = 1 and 79400 
objects for N = 5, k = 2 remain. 

Second catalogue filtering. The raw background 
galaxy catalogues are selected from objects with 
SExtractor isophotal magnitude / > 18.5 and a half-light 
radius which is larger than that measured for stars. The 
resulting catalogues contain altogether ~ 72100 objects 
with iV = 3, fc = 1 and ~ 52900 objects with N = 5, k = 2. 
The principle of PSF correction. The shape of galax- 



ies is influenced by the anisotropic PSF. In order to obtain 
a correct estimate of t he shear 7 from the observed ellip- 
ticity of galaxies e obs , iKaiser et al.l l|l995h developed the 
so-called KSB algorithm. The algorithm relates the ob- 
served ellipticities e obs to the sheared source-ellipticities 
and provides an unbiased estimator of the shear. The cor- 
rection is calculated on the second brightness moments 
Qij of a galaxy with surface brightness 1(6). The quan- 
tity Qij is defined by 

Qij = J d 2 0(9 l -0 l )(0 ] -0 ] )i(e)w (\e-e\ 2 ^ , (31) 

where W is a window function with a smoothing scale 
r g (connected to the object size), which suppresses the 
photon noise of the objects profile at large radii, and 8 
is the centre of the surface brightness. The ellipticity is 
defined as 

Qn - Q22 + 2iQi 2 
e := n Tn • (32) 

Assuming that the intrinsic orientation of galaxies is ran- 
dom, the relation between 7 and e obs in the weak lensing 
regime reads 

7 = (P s )- 1 (e obs -P s Y), (33) 

where P g is the pre-seeing shear polarisabihty which de- 
pends on the smear and shear polarisabihty tensors P sm 
and P sh , and the stellar smear and shear polarisabihty 
tensors p sm * and P sh * . The tensor P s distorts the galaxy 
ellipticity to its true value. Both quantities are calculated 
by means of the observable Qij. The quantity q* is the 
stellar ellipticity (due to the PSF-anisotropy) and is calcu- 
lated from the raw stellar ellipticity e*, q* — (P sm *)~ 1 e* . 
PSF-anisotropy correction. The stars which are used 
for the PSF-anisotropy correction are selected by plotting 
magnitude against half-light radius (mag — plot). All 
stars have the same half-light radius and therefore show 
up as a vertical branch in this plot. Stars which have a 
magnitude of mag = 0.5 lower than the saturated stars 
and which are well above the crowded faint magnitude 
area which contains a mixture of faint stars, galaxies and 
noise detections are selected. Using this sample of stars, 
a third-order two-dimensional polynomial fit with 3.5<r- 
clipping of the stellar ellipticities e* (intrinsically round 
objects are distorted due to the anisotropic part of the 
PSF and are thus elliptical) as a function of position is 
performed. With this, the quantity q* — (p sm *)-! e * a t 
the position of the galaxies is calculated. Fig. 1171 displays 
the stellar ellipticities for all VLT fields before and after 
the PSF anisotropy correction. 

Calculation of P s . The diagonal elements of the P s 
tensor are dominant by a factor of 10 compared to the off- 
diagonal elements and they are approximately equal, so 
that we can estima te P g by P s g l, with P g = 0.5trace[P g ] 
llErben et al.ll200l[) . The stellar smear and shear polaris- 
abihty tensors P sm * and P sh * are calculated for different 
smoothing scales r s . As P s depends on these quantities 
we calculate P 6 according to the galaxy size. 
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Fig. 17. Stellar ellipticities for all VLT fields before (left) 
and after (right) PSF anisotropy correction. 



Third catalogue filtering. All objects having an el- 
lipticity (after PSF correction) of more than 0.8 are re- 
jected. The final catalogues consist of 56800 galaxies for 
N = 3, k = 1 and 42200 galaxies for N = 5, k = 2, re- 
sulting in an average number density of n = 26.3/arcmin 2 
and n = 19.5/arcmin 2 , respectively. The galaxy ellipticity 
dispersion of the final catalogues is on average a e — 0.4. 
In Table seeing condition and galaxy number density of 
the final catalogues are listed for all 50 VLT fields. 
Weighting. Since the corrected galaxy ellipticities are 
very noisy, a weighting scheme according to the noise level 
is introduced. For each galaxy the next twelve neighbours 
are identified in the mag — r^ plane and the variance a\ of 
the ellipticity di s tribut ion of the sub-sample is calculated, 
see lErben et all l)200lh . The variance <r 2 gives an indica- 
tion of the noise level of these galaxies. According to cr 2 
we then determine the weighting factor w as, 

w = (34) 
5.3. Number density and seeing 

The signal-to-noise ratio of a halo detection using the 
aperture mass statistics is proportional to the square root 
of the number density of background galaxies, and the 
expected number density of M ap -peaks strongly depends 
upon the snr. Seeing conditions are therefore crucial. For 
this work we have a large data set of VLT fields, all taken 
under different seeing conditions, so we briefly present 
the dependence of galaxy number density on seeing and 
SExtractor parameter settings. 

Fig. 1181 illustrates the dependence of the number den- 
sity of galaxies used for the weak lensing analysis (fi- 
nal lensing catalogue) on seeing and SExtractor parame- 
ter settings (data from Table 0). The ratio of the num- 
ber density between the SExtractor parameter settings 
N = 3, k = 1 and N = 5, k = 2 is approximately constant 
(n3/n5 « 1.3) as a function of seeing. Using the parame- 
ter settings N = 3, k = 1 we expect the snr of a cluster 
detection to be 1.14 times larger compared to the parame- 
ter settings N = 5, k = 2. However, this is only true if the 
additional sources are not dominated by noise detections 
which do not contain shear information. 




0.4 0.6 0.8 1 ' 0.4 0.6 0.8 1 



seeing ["] seeing ["] 

Fig. 18. The diagram on the left displays the number den- 
sity of galaxies per square arcminutes and its dependence 
on seeing in arcseconds for two different SExtractor pa- 
rameter settings. Solid line: fit to the binned data for 
the parameter settings: N — 5, k = 2; dashed line: 
N = 3, k = 1. The diagram on the right displays the ratio 
between the two number densities caused by the two dif- 
ferent SExtractor parameter settings. The solid line is the 
ratio between the two fits. The number density of galaxies 
are from the final catalogue. Data from Table 

The ratio of the number density of galaxies between 
a seeing condition of 0'.'6 and 0'.'9 is approximately 2.7, 
which corresponds to a ratio of the snr-values of clus- 
ter detections of 1.6. Hence the seeing conditions have a 
tremendous influence on the expected number density of 
detectable clusters. 

We have seen that the galaxy number density strongly 
depends on the seeing conditions and the SExtractor pa- 
rameter settings, and the information about the number 
density of galaxies only makes sense if simultaneously in- 
formation about source extraction and seeing are given. 



5.4. Comparison to numerical simulations 

In this section the results of the simulations are compared 
with real data. Therefore, we apply the M ap -statistics to 
the VLT-images as described in Sect. 12.41 A regular grid 
with a mesh size of 10" x 10" is placed over the data 
fields and the snr is calculated for every grid point. This 
is done for the SExtractor parameter settings 3 contigu- 
ous pixels la above the sky background and two differ- 
ent halo-filter radii, do = 3'.8 and 9o — 5'.7. The number 
density of background galaxies in the numerical simula- 
tions is n = 19arcmin -2 and the ellipticity distribution 
is <j e = 0.34 compared to n < 26 arcmin~ 2 (depending on 
the redshift of the possible galaxy cluster) and o e = 0.4 
for the VLT data. The snr of a halo is proportional to 
V^/cr e . As 719/0.34 = 12.82 « 12.74 = V26/0.4 it is le- 
gitimate to compare the expected number densities of the 
numerical simulations with those of the VLT data. 

In Fig. 1191 we add up all peaks in the snr-maps of 
all VLT fields for the two filter radii, normalise it to one 
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Fig. 19. Comparison of number density between simula- 
tions, including border effects and image treatment (dot- 
dashed lines), and VLT data (solid lines). Dotted lines: 
noise peaks due to randomisation of background galaxies 
of all VLT fields. 

square degree and compare it with the number density of 
peaks obtained in the same way from the 300 subfields cre- 
ated from the numerical simulations (Sect.^). The average 
number of noise peaks per field above a given threshold is 
determined by randomising the orientations of the back- 
ground galaxies 20 times for each field fFig. I19f) . For the 
small filter radius (6q = 3'.8) the number density of peaks 
in the VLT fields is (within the errors) in good agreement 
with the simulations and is significantly above the number 
density of noise peaks. For the large filter radius (6q = 5'.7) 
this is not the case. The number density of peaks in the 
simulations is significantly higher compared to the VLT 
fields and is comparable with that of the noise peaks. The 
reason for this could be that more massive clusters (which 
are better matched by a large filter radius) are present in 
the simulations (as indicated in Sect. l3.3)l than in the VLT 
fields. 

6. Cluster candidates 

In this section the VLT fields are analysed in detail. We 
provide some criteria with which it may be possible to 
separate peaks in the Af ap -map resulting from real halos 
from those peaks resulting from a chance alignment of 
background galaxies. The M ap -statistics are applied to the 
VLT- images using a grid with a mesh size of 3" x 3". We 
use both the polynomial filter (JSJ and the halo-filter (|5J 
with various filter radii (9q S [l',6']) and two SExtractor 
parameter settings. Thereby 30 weak lensing maps (Map- 
maps) are obtained for each of the 50 VLT fields. The 
weak lensing maps are compared with light distribution 
maps. 

The light distribution is calculated on a regular grid, 
where each grid point contains the total flux within a 
weighted circular aperture. The mean light and the stan- 
dard deviation a are calculated from all weighted total 
flux values of the grid points. We optimised the effective- 
ness of the aperture mass filter function, halo-filter U, to 



detected cluster-sized dark matter halos. Assuming that 
light follows the dark matter distribution, we opt for the 
weighting function of the halo-filter U; obtained via Eq. 
(@J and Eq. JJJJ. To calculate the light distribution we use 
galaxies in the magnitude interval I <E [16; 22] and choose 
the same aperture radius as for the weak lensing analysis. 

If we detect peaks in the weak lensing maps with a 
signal-to-noise ratio larger than three which coincides with 
a light overdensity of 2a above the mean light or more, we 
consider these peaks as cluster candidates and perform a 
more detailed analysis of these candidates. We also per- 
form a further analysis if the snr of the weak lensing anal- 
ysis is larger than four. 

In the 50 VLT fields 12 cluster candidates fulfil these 
conditions (in Fig. |23 and Fig. |2S1 all 12 images of the 
candidates with M ap - and light distribution contours are 
presented). In the following, some criteria to substantiate 
or weaken the presumption of the 12 cluster candidates 
to be real clusters are presented. We elaborate only on 
the analysis done with the optimised halo-filter, because 
for the relatively small VLT fields, the polynomial filter 
function is strongly influenced by the border effects (the 
polynomial filter function places weight at large radii; the 
maximum is at 0.7 x 8 , see Fig. 0. For each field the 
number of peaks with snr > 3 for a filter radius of 8q = 3'.8 
is shown in Table El 



6.1. snr-radius plot 

The value of the snr-maximum of the cluster candidates 
is plotted as a function of filter radius, 9q, for different 
filter functions, see Fig. [3U| For this we neglect the fact 
that the position of the maxima can vary by a few pixels. 
A criterion for a promising cluster candidate is that the 
snr should be larger than 3, independent of the halo-filter 
radius in use. Especially for larger radii (in the case of the 
small fields at hand a large radius means 9o ~ 3' — 5') the 
snr should at least be larger than 3. Due to this criterion 
we reject 'vlt36' and 'vlt60' as cluster candidates. 

6.2. Different SExtractor parameter settings 

The 50 VLT fields are analysed for two different 
SExtractor parameter settings, n — 3 contiguous pixel 
k = la above the sky background and a more conserva- 
tive one, n — 5, k — 2. With the first parameter settings 
on the one hand many noisy sources are extracted (many 
of them excluded by the condition that the half-light ra- 
dius of sources must be larger than that of stars), which 
could lower the lensing signal. With the conservative set- 
tings on the other hand, background sources which could 
contribute to the lensing signal are missed. One criterion 
for a promising cluster candidate is that the snr of the 
candidate should be larger than 3 over a large range in 
radius, independent of the SExtractor parameter settings. 
In the case of candidate 'vlt77' (see Fig. I20|l this crite- 
rion is not fulfilled. The snr-curve obtained by using the 
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conservative settings is always lower than 2.6 and is la 
or more below the snr-curve determined with the other 
settings (independent of the filter radius). 



6.3. Exclusion of high ellipticities 

A high snr-peak in the A/ ap -map can be caused by a 
chance alignment of only a small number of galaxies with 
a high tangential cllipticity, especially if the number den- 
sity of background sources is low. Such a statistical fluke 
can be exposed by restricting the absolute value of the 
PSF-corrected ellipticity to |e| < 0.5. We analyse the clus- 
ter candidates once again using the ellipticity restriction 
for the SExtractor parameter settings n = 3, k = 1 (Fig. 
120(1 . The most remarkable example for such a fluke is can- 
didate 'vlt79'. The snr is larger than 4 over a large range 
of radius, independent of the SExtractor parameter set- 
tings. If the galaxy ellipticity is restricted to |e| < 0.5 the 
signal drops below 2a and is therefore rejected. A further 
example is 'vlt42' which is rejected, too. 



6.4. Tangential ellipticity-radius plot 

Another test of the cluster candidates is the tangen- 
tial ellipticity-radius plot. The average tangential shear 
around the snr-maximum of a shear-selected cluster can- 
didate is calculated in rings and is plotted against distance 
from the snr-maximum, see Fig. 1211 Assuming a relaxed 
cluster (an SIS for example), the tangential shear profile 
should follow roughly a 1/ ^-relation. If there are only one 
or two rings which cause the shear signal in the M ap -map, 
then a cluster candidate is rejected. An example is the 
candidate 'vlt54'. 



6.5. Five promising cluster candidates 

We performed an analysis of 50 VLT fields using the halo- 
filter function and obtained weak lensing maps of the 
fields. In 12 of these maps the presence of significant Map- 
peaks is revealed that are associated with overdensities 
in the light distribution (Fig. |H and Fig. Five of 
these remain after a careful analysis (we calculated the 
light distribution, used different source catalogues, anal- 
ysed the fields with various filter radii, excluded high 
ellipticities and calculated the tangential shear profile). 
In Fig. El 3' x 3' clips of the five most promising clus- 
ter candidates (vlt44, vlt86, vlt29, vlt45, vlt55) are dis- 
pl ayed. One of these ca ndidates has been discovered before 
bv iMaoli et"all l|200ll) (candidate vlt55); see ESO Press 
Release 24/00. From the simulations we expect to find 
about two clusters with a snr larger than four using a 
fixed filter radius of 0q = 3'.8. This is comparable with the 
analysed data, as one candidate (vlt55) has a snr larger 
than four using the halo-filter with 9q = 3'.8 and the snr 
of two further candidates (vlt29 and vlt45) reaches four 
for a slightly different filter radius. 



The current data at hand (one colour; /-band; small 
ficld-of-view) allowed us to detect these mass concentra- 
tions, but do not permit a more detailed analysis at this 
stage. As our candidates are not in the centre of the 
FORSl-fields, the gravitational shear and hence the mass 
distribution of the clusters cannot reliably be mapped over 
a large range of radii. With one colour, we cannot obtain 
an estimate for the cluster redshifts or for the redshift 
distribution of the galaxy population used for the analy- 
sis. This, besides an elimination of the contribution from 
foreground galaxies, is essential to obtain accurate mass 
estimates. Hence, to make further progress with these can- 
didates we require multi-colour observations and a larger 
field-of-view. 

7. Summary and conclusions 

We have analysed the ability of the M ap -statistics to de- 
tect massive mass concentrations. We first maximised the 
effectiveness of the M ap -statistics using analytic descrip- 
tions and then applied M ap to synthetic images created 
from numerical simulations. We investigated the influence 
of image treatment and border effects on the snr and num- 
ber density of halos. Finally, we applied M ap to a data set 
obtained with the VLT and compared the results with the 
predictions obtained from the simulations and performed 
a detailed analysis of the cluster candidates. Our major 
findings are as follows. 

We created twelve synthetic images, each covering a 
30' x 30' area, from A-body simulations. We compared 
the halo-filter and the polynomial filter and found that 
4.5 times more peaks with a snr larger than four are de- 
tected with the halo-filter compared to the polynomial 
filter. However, if we take into account the contamination 
ratio of noise peaks to the total number of peaks in the 
Af ap -maps, the difference in efficiency to find real clusters 
for the two filter types is much less distinct. 

We studied the effect of weighting and image treatment 
on the snr of peaks in the weak lensing map. We found 
that, on the one hand, the image treatment lowers the snr 
significantly so that the expected number density of halos 
decreases by a factor of two. On the other hand, weighting 
has only a weak influence on the snr. To compare real 
data with simulations, both effects have to be taken into 
account. Based on our findings in Sect. l3~5l we concluded 
that it is more efficient for future weak lensing surveys to 
propose for medium deep images than for a few very deep 
images. 

We studied the influence of the border effects on the 
number density of peaks in the weak lensing maps. For 
that, we subdivided the fields of the numerical simulations 
into 300 subfields (all having the same size as the VLT 
fields). We found that the border effects affect the snr of 
clusters significantly and consequently has a large impact 
on the number density of peaks in the M ap -maps. The 
number density of peaks in the M ap -maps decreases by 
a factor of two. Including the border effects, we expect 
for VLT-sized images ~ 3 peaks per square degree with 
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a snr larger than four using the halo filter with a filter 
radius of 8 n = 3'.7, see Table In the future, ground- 
based wide-field images will be common, so this effect will 
not play an important role. However, for future space- 
based missions with smaller field of views (compared to 
ground-based cameras) this effect still has to be taken into 
account. 

We then performed a statistical peak analysis of the 50 
VLT fields and found that the number density of peaks (for 
snr > 3 and a filter scale of 3'.8) is comparable with the 
number density obtained from the numerical simulations, 
taking into account border effects and image treatment. 
However, for a filter scale of 5'.7 the number density of 
peaks in the VLT fields is significantly smaller compared 
to the simulations and is comparable with the number 
density of noise peaks, see Fig. We point out again 
that the reason for this could be that low-mass clusters 
are present in the VLT fields which are not matched with 
the larger filter function. If this would be the case, then 
this is a way to constrain the dark matter halo size in a 
statistical way. 

Finally, we reported the results of a detailed analysis 
of the VLT fields. We detected several M ap -peaks with 
a snr > 3.0 in our 50 fields. For 12 of the M ap -peaks 
we could associate an overdensity in the light distribu- 
tion. These cluster candidates were analysed in detail by 
using different filter radii, SExtractor parameter settings 
and ellipticity cuts. Furthermore, the tangential ellipticity 
and its dependence on distance to the snr-maximum was 
analysed. Finally, five promising candidates remain after 
selection which need a follow-up observation in different 
filters to clarify their nature. 

In this work we are quite sceptical on the efficiency 
of finding individual clusters with masses less then M w 
3 x 10 14 Mq for redshifts larger than z = 0.3. But since 
there is much more information in low-mass clusters than 
in high mass ones, simply because they are so much out- 
numbered, the goal should not be a blind search for in- 
dividual clusters, but a statistical blind search, especially 
for future large weak lensing surveys. A statistical peak 
analysis has already b een applied to observable data by 
iMivazaki et al. | ll2002h and is a valuable tool to explore 
cosmological models. The major part of our work should 
be seen as a foundation stone for the comparison of weak 
lensing surveys with ray-tracing through V-body simu- 
lations since we now consider image treatment or other 
observational effects (like the border effect). This will 
improve the manner by which cosmological models can 
be explored with galaxy clusters. With better synthetic 
data (more realistic background galaxies could be sim- 
ulated with shapelets, for instance) generated from ray- 
tracing through A^-body simulations for different cosmo- 
logical models and a much larger survey (compared to 
our 0.64 square degree) it would then be possible to con- 
strain cosmological parameters and especially the statis- 
tics about the distribution of dark matter. 



Appendix A: Calculation of the number density of 
halos 

We first describe in this appendix in detail the calcula- 
tions to obtain the aperture mass, M ap , given a trun- 
cated NFW-profile, a redshift distribution of background 
galaxies and the new filter function (halo- filter). We show 
then how we calculate the expected number density of ha- 
los assuming a halo distribution which utilise the Press- 
Schechter theory. The calculations are restricted to a flat 
universe (fi m + f^A = 1)- 

The aperture mass, M ap , is given by 



M a 



dz s k(&, z\, z s ) U {&) p(z s ), 



(A.l) 



where U is the halo- filter obtained by using Eq.^Jand Eq. 
|8J p is the redshift distribution of galaxies given in Eq. fH] 
The convergence, k, is that of an N FW-profile and 
is trun cated at the virial radius, r v i r , see iTakada fc .Tainl 
( 2003) and is given by 



n(y) = 3tt m 5 s r s 



Ho 
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(A.2) 



with y = r/r s , where r s is the scale radius, D denotes the 
angular diameter distance and the function / is given by, 



f(y) 



where 



A + (1 - y 2 )' 3 / 2 arccosh^T±^ 

Cw-l 2+CN. 



A = 



y 



3(l+c N ) l+c N ' 

A - (y 2 - 1)~ 3 / 2 arccos 
0; 



c N 



»(1+cn) ' 

j) 2 +CN . 

2/(l+c N ) ' 



y < 1 

1 < y < c N 
y > cn 



2/(1 + c N ) 
The quantity S s is given by 
A . 

5 S = 



(A.3) 



(A.4) 



3 log(l + c N ) - c N /(l + c N ) ' 
The concentration parameter of an NFW profile is 

cn = — (A.5) 

r s 

and can also be expressed by ijBullock et al.l F2001) 
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where we set c* = 8 for an open universe. The mass within 
a sphere of radius r v ; r (virial radius) is 



log(l + cn) 



c n 



1 



CN 



(A.7) 



The virial mass can also be defined by the spherical top- 
hat collapse model as 



-in 



S viY (z)p rl i 



3 

vir ' 



(A.6 



18 



M. Hetterscheidt et al.: Searching for galaxy clusters using the aperture mass statistics in 50 VLT fields 



with 



PO — Pcrit^mj 
i Pr.rit 

reads llBullock et al.ll20Q]|) 



(A.9) 



where p rr \t is the critical density. The virial overdensity 



5 vil (z) = (18tt 2 + 82a: - 39a; 2 ) 



where x = f2(z) — 1 and 



(1 + Zi) 3 fi m + (1 + Zl) 2 (l - fim - fi A ) + 



(A.10) 
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To calculate the expected number density of halos per 
steradian with ape rture mass larger than AL n we use the 
formula derived bv lKruse fc Schneiderl |^}99), 



A^(> M a p) 

with 

G{z h M ap ) 
and 



n l ^i^7^Df( Zl )G(z h M ap ), (A.12) 
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and 



N halo (M, z x ) dMdV c = J\ -L ^ila 



M <t 2 (M) 



der(M) 
dM 



x exp(-g^g}d:i/<ii; 



(A.15) 



is number of objects in the comoving volume dV c with 
mass in the interval dM . The quantity 



5 C Tit(z\) = 



u crit 



D+(zi,£1 ,£Ia) 



(A.16) 



is the critical density threshold for spherical collapse which 
depen ds on the linear growth factor, D + ijLacev fc Cold 
Il993l) . The quantity cr(M) is the present linear theory rms 
density fluct uation with shape paramet er T and normali- 
sation <7 8 . As lKruse fc Schneiderl l|l999|) we use the fitting 
formu lae A6-A19 given in the appendix of iNavarro et alJ 
l|l997(l to calculate u(M) and <5°. it . 
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Fig. 20. The detection significance of the shear-selected cluster candidates. Shown is the maximum signal-to-noise 
ratio of the cluster candidates in the M ap -map and its dependence on filter radius 9o for different filter (halo- and 
polynomial filter function) and SExtractor parameter settings (n = 3 and n = 5 contiguous pixel, k = la and k = 2a 
above the sky background). To exclude that the snr of a cluster candidate is dominated only by a small fraction of 
galaxies with a high ellipticity we also show a M ap -analysis for which we exclude ellipticities with |e| > 0.5. The capital 
X denotes the promising cluster candidates. 
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Fig. 21. The average tangential shear calculated in rings and its dependence on the distance, 9 (in arcmin), from the 
snr-maximum of the shear-selected cluster candidates. The error bars are calculated from the cross component of the 
shear. The solid squares without error bars are the cross components of the shear and are shown as a comparison. 
The numbers in the upper part of each figure denote the number of galaxies per ring. Note that some of the cluster 
candidates are close to the edge of a field. 
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Fig. 22. VLT fields of all cluster candidates (indicated by circles). The M ap -contours are black (snr-values are 1.5, 
2.0, 2.5, 3.0, 3.5, 4.0) and the light distribution contours are white (snr-values are 1, 1.5, 2, 2.5, 3, 3.5, 4). The capital 
X denotes fields with promising cluster candidates. 
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Fig. 23. As Fig. [521 for the other six fields with cluster candidates. 



22 M. Hetterscheidt et al.: Searching for galaxy clusters using the aperture mass statistics in 50 VLT fields 




22*30^42* 40* 3fl* 3fi* S4* 32* £2^1 5° IS 3 12 1 »' 0&* 



Fig. 24. 3' x 3' clip of the most promising cluster candidates. The M ap -contours are black (snr-values are 1.5, 2.0, 
2.5, 3.0, 3.5, 4.0) and the light distribution contours are white (snr-values are 1, 2, 3, 4, 5). 
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Table 5. List of all 50 VLT fields. n5 and n3 are the number density of galaxies per square arcminute in each field of 
the final catalogue with SExtractor parameter settings N — 5 and N — 3 connected pixels with k = 2a and k = la 
above the mean background value, respectively. The seeing given in this table correspond to the stacked image. iV(> 3) 
denotes the number of peaks per field with a larger snr than 3. For the cluster candidates the maximum snr and the 
corresponding halo-filter radius are listed. The maximum snr of the rejected candidates are in brackets. 





RA ( , T2000) 


DEC r T20001 




Ti5 


7l3 


Tad [mfl cl 


N(> 31 


' 11 max V^O / 


vlt27 


nn 59 28 i 

ww o o wu. i 


—00 18 28 

WW -LO ZjO 


0'.'72 


8 


15 


24.5 


3 




vlt28 


ni 31 4n 3 

W -L O-L 4:W.O 


— nn 22 28 

WW ZiZi ZiO 


q"54 


27 


40 


24.9 


o 






01 df\ R 

U-L Of L ±yJ.O 


— nn n^ ^i 

UU UO Ol 


n"4Q 

u. 4y 


^9 


OU 


25.1 


2 


4 4 f t i- 

4:. 4: ^O.O ^ 


vlt30 


09 98 44 


— nn n^ 9<s 

UU UO ^u 


U. 04: 


9^ 
^o 


OO 


24.9 


Q 
O 




vlt31 


m nn 9i 8 

U _L uu z< _i_ . o 


— n3 1 5 31 

WO -LU O-L 


q"57 


24 


35 


25.0 


5 




vlt33 


n9 nn 08 i 

IJU UO . -L 


— n3 nn 3n 

WO WW ow 


q"44 


29 


44 


24.9 


5 




vlt35 


nn 59 35 3 

uu Otj oo.o 


— n6 in n5 

WW 1U wo 


0"73 


17 


24 


24.8 


3 




vlt36 


91 98 53 1 

U _L ZiO UO . -L 


— n6 ni 39 

WW U-L Ot? 


0'.'68 


25 


26 


24.8 


3 


f4 3 Cl'4n 


vlt37 


ni 57 n5 8 

u _L u J uu . o 


—96 n5 m 

WW WO W-L 


0'.'90 


14 


18 


24.5 


I 




vlt39 


91 3n 45 3 

Zi_L OU 4:U.U 


— 99 58 45 

VJiJ OO 4:0 


0'.'76 


18 


23 


24.8 


I 




vlt40 


99 n4 37 9 

£t£i U4: O 1 . ij 


—in 15 n9 

_LU 1U Uf 


0'.'71 


19 


27 


24.9 


3 




vlt42 


99 9Q 9Q 9 


—in 12 m 


n"79 

u. / ^ 


1 s 


1 7 


94 7 


2 


[o.o (.o.ej] 


vlt43 


91 »n 9. » 


— 15 11 48 


n"55 

u. OO 


9n 

^u 


9Q 


9<t n 

^o.u 


4 




vlt44 


99 09 16 6 

WZ -LW.W 


— 14 53 n3 

J-rfc OO WO 


0'.'64 


24 


25 


24.9 


2 


4.4 (l'.4) 


vlt45 


22 3n 41 8 


— 14 54 55 


0'.'46 


28 


32 


24.8 


5 


4 1 (^l 


vlt46 


22 m 42 2 


— 2n in 55 

ZW -LU OO 


0'.'65 


24 


33 


25.0 


4 




vlt47 


22 29 33 8 


— 9n 14 44 

iiW It TtTt 


0'.'51 


25 


38 


25.0 






vlt48 


91 ^n 4^ ^ 

£L OU 4:0.0 


— 94 , , 4n 

Zf4: OO 4:W 


U. UO 


22 


nn 

OU 


24.8 






vlt49 


21 58 44 7 

Z _L OO tt . I 


—24 57 15 


0'.'62 


21 


30 


24.6 


3 




vlt50 


22 3n 43 8 

OW 4:0.0 


—25 ni 42 

ZiO W-L Tilj 


0'.'55 


24 


36 


25.0 


I 




vlt51 


2n 59 39 5 

OiJ OW.U 


— 3n 18 31 

OW -LO O-L 


0'.'62 


15 


22 


24.5 


I 




vlt52 


99 nn 9(S 9 


— »n ni 4 . 

OU U-L 4:0 


U. UO 


1 Q 


9fi 
^u 


24.8 


I 




vlt53 


22 31 15 3 

O _L _LO.O 


— »n n7 i , 

OU U 1 -LO 


0'.'55 


27 


39 


24.7 


o 




vlt54 


21 29 53 6 

Z; _1_ ZjU OO.O 


— 34 51 52 

04: O -L 


0'.'57 


22 


24 


25.1 


5 


f3 7 C^'M 

[O. / (4 )\ 


vlt55 


99 nn 14 i 

WW It:. -L 


— 35 3n 54 

OO OU U4: 


q"53 


22 


33 


24.6 


4 


4.3 (3'.8) 


vlt56 


22 3n nf! 4 

OU WW.rr 


— 35 in 33 

OO -LU OO 


0'.'83 


13 


19 


24.5 


I 




vlt57 


91 98 f\d Q 


_QQ AO n9 

oy 4:y uzr 


n"55 

U. OO 


24 


OO 


9^ n 

^o.u 


Q 
O 




vlt58 


22 nn nfi 7 

ZiZi WW WW . i 


— 4n n4 , , 

4:U Ut: OO 


q"49 


26 


37 


24.8 






vlt59 


22 29 1 1 8 

ZjZj ZiU _L _L . O 


— 39 36 28 
oy ou 


0'.'70 


21 


27 


24.7 






vlt60 


22 59 24 4 

Oi7 


—in ni 29 

_L W W-L ZiU 


q"47 


28 


33 


25.0 


4 


\A 1 Cl'411 

[4.1 


vlt61 


22 59 24 2 


— 15 08 47 


0"47 


34 


47 


25.0 






vlt62 


22 59 91 8 

ZiZi kj U W -L . O 


— 19 44 03 


q"47 


28 


40 


25.0 






vlt63 


22 59 39 5 


—24 52 51 

O-L 


q"49 


23 


35 


25.0 


5 




vlt64 


22 59 56 1 

ZiZi KJ U OW.-L 


—30 14 27 


0'.'60 


19 


27 


24.5 


2 




vlt65 


23 nn 44 3 

ZO WW 4:4:. O 


— 34 55 26 

Ot: OO 


q"54 


25 


38 


25.0 


4 






9 Q ni 94 & 

ZrO U-L Zr4:.0 


—40 9 . 9n 

4:U ZrO ZrU 


n"77 

U. / 1 


12 


1 7 


9^ fi 

^o.u 


2 




vlt75 


21 28 14.7 


-20 07 18 


0'.'56 


26 


36 


24.9 


4 




vlt76 


21 32 21.0 


-30 25 57 


0'.'63 


19 


26 


24.6 


3 




vlt77 


14 59 07.4 


00 07 54 


fy.'80 


16 


16 


24.8 


2 


[3.6 (3'.4)] 


vlt78 


14 59 03.2 


05 11 32 


rt'so 


31 


34 


25.0 


1 


[3.2 (5')] 


vlt79 


14 59 32.7 


10 13 19 




13 


15 


25.0 


3 


[4.5 (5')] 


vlt80 


15 30 17.5 


00 10 58 


fy.'60 


23 


32 


24.8 


3 




vlt81 


15 29 40.4 


04 54 10 


0763 


20 


27 


24.7 


1 




vlt82 


15 28 59.7 


10 14 59 


0'.'59 


22 


31 


24.6 


4 




vlt83 


15 59 00.7 


-00 07 14 


07 87 


11 


16 


24.3 


4 




vlt84 


16 03 35.0 


05 10 46 


(/.'91 


11 


15 


23.8 







vlt85 


15 56 47.6 


10 17 28 


0766 


19 


27 


24.7 







vlt86 


16 00 30.1 


14 58 35 


0778 


16 


16 


24.6 


1 


4.0 (1'.7) 



24 



M. Hetterscheidt et al.: Searching for galaxy clusters using the aperture mass statistics in 50 VLT fields 



References 

Bartelmann, M. & Schneider, P. 2001, Phys. Rep., 340, 
291 

Bertin, E. & Arnouts, S. 1996, A&AS, 117, 393 
Braincrd, T., Blandford, R. D., & Smail, I. 1996, ApJ, 
466, 623 

Bullock, J., Kolatt, T., Sigad, Y., ct al. 2001, MNRAS, 
321, 559 

Dahle, H., Pedersen, K., Lilje, P. B., Maddox, S. J., & 

Kaiser, N. 2003, ApJ, 591, 662 
Dietrich, J. P., Schneider, P., Clowe, D., et al. 2004, 

preprint astro-ph/0406541 
Erben, T., van Waerbeke, L., Bertin, E., Mellier, Y., & 

Schneider, P. 2001, A&A, 366, 717 
Erben, T., van Waerbeke, L., Mellier, Y., et al. 2000, A&A, 

355, 23 

Hamana, T., Takada, M., & Yoshida, N. 2004, MNRAS, 
350, 893 

Hcnnawi, J. F. & Spergel, D. N. 2005, preprint astro- 
ph/0404349 

Jain, B. & van Waerbeke, L. 2000, ApJ, 530, LI 
Jenkins, A., Frenk, C. S., White, S. D. M., et al. 2001, 

MNRAS, 321, 372 
Kaiser, N., Squires, G., & Broadhurst, T. 1995, ApJ, 449, 

460 

Kaiser, N., Squires, G., Fahlman, G., & Woods, D. 1994, in 

Clusters of galaxies, proc. XlVth Moriond astrophysics 

meeting, Mcribcl, France, 1994, 269 
Kruse, G. & Schneider, P. 1999, MNRAS, 302, 821 
Lacey, L. & Cole, S. 1993, MNRAS, 262, 627 
Maoli, R., van Waerbecke, L., Mellier, Y., et al. 2001, 

A&A, 368, 766 
Maturi, M., Meneghetti, M., Bartelmann, M., Dolag, K., 

& Moscardini, L. 2004, preprint astro-ph/0412604 
Miyazaki, S., Hamana, T., Shimasaku, K., et al. 2002, 

ApJ, 580, L97 
Navarro, J., Frenk, C, & White, S. 1996, ApJ, 462, 563 
— . 1997, ApJ, 486, 493 

Press, W. & Schechter, P. 1974, ApJ, 187, 425 
Schirmer, M. 2004, PhD thesis, Universitat Bonn 
Schirmer, M., Erben, T., Schneider, P., ct al. 2003, A&A, 

407, 869 
— . 2004, A&A, 420, 75 
Schneider, P. 1996, MNRAS, 283, 837 
Schneider, P., Van Waerbeke, L., Jain, B., & Kruse, G. 

1998, MNRAS, 296, 873 
Takada, M. & Jain, B. 2003, MNRAS, 340, 580 
Umetsu, K. & Futamasc, T. 2000, ApJ Lett., 539, L5 
Van Waerbeke, L. 2000, MNRAS, 313, 524 
Von der Linden, A., Erben, T., Schneider, P., & 

Castander, F. J. 2005, preprint astro-ph/0501442 
White, M., Van Waerbeke, L., & Mackey, J. 2002, ApJ, 

575, 640 

Wilson, G., Kaiser, N., Luppino, G. A., et al. 2001, ApJ, 
555, 572 

Wittman, D., Margonincr, V. E., Tyson, J. A., et al. 2003, 
ApJ, 597, 218 



Wittman, D., Tyson, J. A., Margoniner, V. E., et al. 2001, 
ApJ, 557, 89 



